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Abstract. We study the properties of an TV-site tight-binding ring with parity and 
time-reversal (PT) symmetric, Hermitian, site-dependent tunneling and a pair of non- 
Hermitian, PT-symmetric, loss and gain impurities ±17. The properties of such 
lattices with open boundary conditions have been intensely explored over the past 
two years. We numerically investigate the PT-symmetric phase in a ring with a 
position-dependent tunneling function t a (k) — [k(N — k)] a / 2 that, in an open lattice, 
leads to a strengthened 'PT-symmetric phase, and study the evolution of the PT- 
symmetric phase from the open chain to a ring. We show that, generally periodic 
boundary conditions weaken the 'PT-symmetric phase, although for experimentally 
relevant lattice sizes N ~ 50, it remains easily accessible. We show that the chirality, 
quantified by the (magnitude of the) average transverse momentum of a wave packet, 
shows a maximum at the PT-symmetric threshold. Our results show that although 
the wavepacket intensity increases monotonically across the PT-breaking threshold, 
the average momentum decays monotonically on both sides of the threshold. 



1. Introduction 

One of the founding tenets of quantum theory is that the Hamiltonian of a "closed" 
system is Hermitian with respect to the standard inner product; this property of the 
Hamiltonian implies that the energy spectrum of the system is purely real, and that 
its eigenvectors form a complete orthonormal set [U 121 E]- The quantum theory, 
based on this and other axioms, has been phenomenally successful in explaining 
experimental results, and predicting various non-trivial and counterintuitive phenomena. 
Traditionally, non-Hermitian Hamiltonians, with eigenvalues that have an imaginary 
part, have been used to model open, dissipative quantum systems such as a resistor [I]. 
Since the seminal work of Bender and co-workers in the late nineties, over the past 
decade, a new class of non-Hermitian Hamiltonians has been extensively explored [HI El 
[7]. Although these Hamiltonians are not Hermitian under the standard inner product, 
they are invariant under the combined parity (V) and time-reversal (T) operation, and 
thus are called PT-symmetric Hamiltonians. The eigenvalues e n of such a Hamiltonian 
are real over a range of parameters, although the corresponding eigenvectors |n) are 
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not orthogonal under the standard inner product. This parameter space where all 
eigenvalues are real is called the PT-symmetric phase, and the emergence of the first pair 
of complex eigenvalues that occurs when the Hamiltonian leaves this parameter space is 
called PT-symmetry breaking. Since time-reversal is an anti-linear operator, it follows 
that when the VT symmetry is broken, the eigenvectors are no longer simultaneous 
eigenfunctions of the combined VT operation. 

Over the past three years, PT-symmetric Hamiltonians for tight-binding lattice 
models have been extensively investigated [HI El HD1 EH]- They have included one 
dimensional chains with non-Hermitian, PT-symmetric tunneling [12], non-Hermitian 
PT-symmetric pair of on-site impurity potentials [13], two-dimensional lattices with 
PT-symmetric potential [14] . one-dimensional chain of dimers with PT-symmetric 
impurities on each dimer [T5], and PT-symmetric spin models [TB] where the effects of 
boundary conditions are claimed to be negligible. These theoretical investigations have 
been accompanied by the experimental exploration of consequences of VT symmetry 
breaking in coupled optical waveguides [TT] HE1 HE] and electrical circuits [20J. These 
ongoing studies have hinted at the deep and rich set of phenomena that follow from non- 
Hermitian Hamiltonians with balanced loss and gain. Although most of these studies 
have focused on open chains, some recent works have shown that the effect of boundary 
conditions on the PT-symmetric phase is non-trivial [2TI |2"2"] . 

The current and emerging promising candidates for PT-symmetric systems - 
coupled optical waveguides and coupled electrical circuits - are presently focused on 
the realization of PT-symmetric open chains. For electrical circuits, since the spatial 
circuit geometry can be easily manipulated, introducing inductive coupling between the 
"end sites of an open chain" is relatively straightforward [20J. On the other front, 
two dimensional (2D) lattices of coupled optical waveguides have been extensively 
experimentally explored [231 121] in the context of Anderson localization. Thus, creating 
a one-dimensional system with periodic boundary conditions - a 2D lattice with only 
boundary waveguides and no interior waveguides - is feasible with current fabrication 
technology. Therefore, realization and experimental investigation of one dimensional 
PT-symmetric rings is likely to occur in near future. 

In this paper, we numerically investigate the PT-symmetric phase diagram in 
a one dimensional system with periodic boundary condition (a PT-symmetric ring). 
We consider an iV-site lattice with a position-dependent, parity-symmetric tunneling 
function t a (k) = t [k(N -k)] a/2 , and a pair of balanced gain and loss impurities ±£7. In 
the Hermitian limit, 7 = 0, such non-uniform open chain shows tunable energy spectra, 
wave packet evolution, and Hanbury-Brown-Twiss correlations [22]. Our primary results 
are a follows: i) For a > 0, we find that the PT-symmetric phase is generally weakened 
in comparison with its open chain counterpart. However, for experimentally relevant 
system sizes N ~ 50, the critical impurity strength jpt is of the same order of magnitude 
as its open-chain counterpart, ii) For a < 0, 7pt for a ring is substantially enhanced from 
its open chain value for all impurity locations, iii) The average transverse momentum of 
a wave packet shows a maximum at the PT-symmetry breaking threshold, and decays 
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monotonically on the two sides of it. 

The plan of the paper is as follows. In Sec. [2j we introduce the tight-binding model 
and discuss the differences between properties of an open chain and a ring that become 
relevant for small N. In Sec. [3j we present the results for the PT-symmetric phase 
diagram as a function of the scale factor < A < 1 that determines the tunneling 
between end points of an open chain. In Sec. |4j we discuss the time evolution of a wave 
packet across the "PT-symmetric threshold. We show that the chirality, which encodes 
the average clockwise or counterclockwise motion of the wave packet around the ring, 
is quantified by the average momentum and shows a maximum at the threshold. The 
paper is concluded with a brief discussion and open questions in Sec. |5j Throughout 
this paper, we consider coupled optical waveguides with complex refractive index and 
gain as the prototype realization of a "PT-symmetric ring [TTJ [TBI HH] ■ 



2. Tight-Binding Model 

We start with a one-dimensional, tight-binding, iV-site lattice with two non-Hermitian 
impurities (+i r y, — ij) located at mirror-symmetric sites (m,rh), where 1 < m < N/2 
and rh = N + 1 — m > m. The Hermitian, position-dependent tunneling Hamiltonian 
for the open chain is 

JV-l 

H = - ^ ta{i){(4+l a i + a l a i+l)> (!) 
i=l 

where t a (k) = t [k(N — k)] a ^ 2 = t a (N — k) > is the tunneling between sites k and 
k + 1, and a\ (a^) represents the creation (annihilation) operator for a single-particle 
state \k) localized on site k = 1, . . . , N. The difference equation for the coefficients of an 
eigenfunction \ip) = X^fcLi/fcl^) °f the Hamiltonian Eq.|l| with eigenenergy e is given 
by 

t a (k)f k+1 + t a (k - l)/ fc _i = -e/ fc , (2) 

where open boundary conditions imply t a (N) = = t a (0). We remind the reader 
that since the tunneling amplitudes is site dependent, traditional methods for solving 
Eq.(j2]) such as the Bethe ansatz are not applicable. Apart from a few exceptions [TU] . 
the eigenvalues and eigenvectors of Eq. ([T]) must be obtained numerically. The impurity 
potential is given by 

V = i^(a ] m a m - a^am) ^ V ] '. (3) 

The action of the parity operator is given by Va k V = at, where k = N + 1 — k and the 
anti-linear time-reversal operator implies TiT = —i. From these definitions, it follows 
that the non-Hermitian potential V is odd under parity and time-reversal individually, 
and but obeys VTVVT = V ^ . The "PT-symmetry breaking in such a non-uniform 
open chain has been extensively studied, albeit only numerically due to the absence of 
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analytical methods that are applicable to a system with site-dependent tunneling [26J. 
The Hamiltonian for a "PT-symmetric ring is given by 

H(X) = H + V + \t R (a\a N + a^aj) , (4) 

where we choose the tunneling between the end points of the open chain, sites 1 and 
N, as tR = t a (l) = t a (N) and the scale-factor < A < 1 allows us to continuously 
extrapolate from an open chain to a PT-symmetric ring. Before discussing the results 
for a ring, A > 0, we briefly recall the results for an open chain and establish the 
terminology [26]. In the PT-symmetric phase, the eigenvalue spectrum of an iV-site 
open chain is non-degenerate and symmetric about zero, and the eigenfunctions for 
energies ±e are related to each other. The bandwidth A' a of the energy spectrum 
scales as A' a ~ N a for a > and A' a ~ jV~H/2| for a < because the bandwidth is 
determined by the largest tunneling amplitude. We use quarter-bandwidth, A a = A'^/4 
as the energy scale and note that when a = 0, the threshold impurity strength is given 
by 7 PT /A a=0 = 1. 

Note that although the distinction between an open chain and a ring is expected 
to vanish [TU] in the limit N — > oo, for small N the differences between the two can be 
substantial. As an extreme case, let us consider a uniform 3-site lattice with nearest- 
neighbor tunneling £ - The non-degenerate, particle-hole symmetric spectrum of such 
an open chain is given by E = { — y/2t , 0, +\/2t }. On the other hand, the spectrum of 
a three-site ring is given by E = {— 2t Q , to, to} and is, in general, asymmetric about zero 
and degenerate. 
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Figure 1. The left-hand panel shows the spectra for Hamiltonian [T] with N = 500 
and a = { — 1,0,1,2}, where the energy is normalized by its maximum value. The 
right-hand panel shows the difference between energy spectra for a ring and an open 
chain for an N — 30 lattice. When a > the spectral differences are most pronounced 
near the center of the band. For a < 0, the spectral differences are greatest at the band 
edge and represent the changes that occur in eigenstates localized at the two ends of 
the open chain |25] . These results show that energy spectrum of a 'PT-symmetric ring 
is different from that of an open chain for experimentally relevant lattice sizes. 
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The left-hand panel in Fig. [T] shows typical spectra for Hamiltonian[T]with iV = 500 
and a = { — 1,0,1, 2}. The energy is normalized by its maximum value and since N ^> 1, 
the spectrum is virtually identical for an open chain and a ring. When a = (black 
line) we obtain the expected cosine spectrum, for a = 1 (green line) and a = 2 (red 
line), the spectrum is quasilinear, and when a = — 1 (blue line), the spectrum band- 
edges represent eigenstates localized at the end of the chain that are generically present 
when a < [23] . The right-hand panel in Fig. [T] shows the difference between energies 
AE = -Ering — -Echam of a ring and an open chain for a lattice with N = 30 sites and 
a = {0, 1, —1}. We find that for a > the difference is greatest near the center of 
the band, whereas for a < the spectral difference is greatest at the band edges. This 
is expected since the localized edge-states that exist for a < are most influenced by 
the introduction of periodic boundary condition. The spectral differences vanish with 
increasing N but remain pertinent for experimentally relevant lattice sizes {271 |28j 129] . 

3. PT-symmetric Phase Diagram 

We now discuss the numerically obtained PT-symmetric phase diagram for the 
Hamiltonian H(X) as a function of increasing loss and gain impurity strength 7, 
fractional location of the gain impurity, li = m/N < 1/2, and the dimensionless scale 
factor A that extrapolates between an open chain (A = 0) and a ring (A = 1). The VT- 
symmetric phase is called robust if the critical impurity strength measured in units of the 
energy scale, ^ PT /A a (N) is nonzero as iV — >• 00; it is called fragile if 7pt/ A a (N) — > 
as N 00. In an open chain li = 1/N corresponds to farthest impurities whereas the 
in a ring, farthest impurities correspond to li ~ 1/4. We remind the reader that in an 
open chain, when a > the critical value of impurity strength 7pt(/u) decreases as the 
distance between the impurities increases whereas for a < 0, the PT-symmetric phase 
is vanishingly small for almost all values of impurity location li [26J. 

Figure [2] shows the typical evolution of the PT-symmetric phase diagram in the 
(j/A a ,n) plane as a function of the scale parameter A. These results are for a lattice 
with N = 30 sites, and a = 1 (left-hand panel) and a = 2 (right-hand panel). For A = 
(black circles), the PT-symmetric phase is robust when the loss and gain impurities 
are closest to each other, li = 1/2. Then the critical impurity strength is given by 
j PT /A a (N) = 1 when N is even and ^ PT /A a (N) = 1/2 when N is odd [26J. For A > 0, 
we see that the critical impurity strength ^p T is, in general, suppressed relative to its 
value for an open chain for all values of impurity positions li. Thus, the PT-symmetric 
phase in a ring with non-uniform tunneling is weaker than its counterpart in an open 
chain. However, results in Fig. [2] also show that the critical impurity strength jpj- is still 
an appreciable fraction of its value in an open chain. The left-hand panel shows that for 
a lattice with a = 1, the critical impurity strength 7pt(/^) is strongly suppressed except 
for a few specific values of /1; this feature is robust irrespective of N and is related to the 
exactly linear spectrum of an a = 1 open chain [TOj [12] . The right-hand panel shows 
that for a = 2 lattice, with a non-linear spectrum, the critical impurity strength reaches 
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a = 1 a =2 
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Figure 2. The left-hand panel shows the "PT-symmetric phase diagram as a function 
of impurity strength 7/A a and fractional impurity position /i = m/N for an ./V = 30 
site lattice with a = 1; for 7 > 7pp, the eigenvalues of the non-Hermitian Hamiltonian 
Eq.Q become complex. The right-hand panel shows corresponding results for the 
a — 2 lattice. In both cases, the PT-symmetric phase is maximally robust at /i = 1/2 
for an open chain, A = (black circles). As A is increased, thus increasing the tunneling 
between sites 1 and N, the critical impurity strength 7pt(m) remains essentially 
unchanged from its open chain value for A < 0.5. As A is increased further (red 
squares, green diamonds), the PT-symmetric phase in the ring (A = 1, blue triangles) 
is weakened for all impurity positions. Thus, a minor change in the Hamiltonian, 
Eq.Q leads to a suppression of the critical impurity strength 7pt(m) even when the 
impurity location is far away from this change. 



a plateau 7pp/A Q ~ 0.3 for most impurity locations /1 > 0.2. As A is increased from zero 
(red squares, green diamonds) to one (blue triangles) the PT-symmetric Hamiltonian, 
Eq.Q, is perturbed only slightly. The tunneling introduced between sites 1 and N, 
\t a (l) = Xt a (N), is the smallest among all links for a > 0. However, it has a dramatic 
effect on ^ PT even when the loss and gain impurities ±27 are away from this tunneling 
link, fi > 0.1. 

When a < 0, the tunneling introduced between sites 1 and N is the largest among 
all links and is proportional to the bandwidth A a (N) ~ j\M a '/ 2 . Therefore, we can 
expect that the fragile "PT-symmetric phase in an open chain will be strengthened 
in a ring. Fig. [3] shows the "PT-symmetric phase diagram for an JV = 30, a = —1 
lattice. When A = (black circles), we obtain the extremely fragile phase diagram of 
an open chain. As the tunneling is increased from A = 0.2 (red squares), A = 0.7 (green 
diamonds) to A = 1.0 (blue triangles), we see that the VT -symmetric phase in a ring is 
substantially strengthened for all impurity positions. Since PT-symmetric rings (A = 1) 
with a > and a < can be mapped onto each other with an appropriate redefinition 
of the impurity location fi, the critical impurity strength 7pt(a*) decreases with /i when 
a < (Fig. [3]) whereas it increases with fx for a > (Fig. [2]). 

We remind the reader that in a ring with constant tunneling, accessible to analytical 
treatment, the critical impurity threshold is zero [22]. These numerical results show 
that in a ring with non-uniform tunneling t a (i) there is a large region, below the blue- 
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Figure 3. PT-symmetric phase diagram for an TV = 30 lattice with a = —1, as a 
function of the scale-factor A that determines the tunneling between sites 1 and N 
and extrapolates from an open chain (A = 0) to a ring (A = 1). When A = (black 
circles) , the 'PT-symmetric phase is fragile everywhere except when the impurities are 
closest, /i = 1/2. For A > 0, the dimensionless critical impurity strength jpt/Ao is 
significantly enhanced for all impurity locations fx, and, in contrast with the a > 
case, 7pt(m) f° r the ring decreases with increasing \i. These results are expected since 
the tunneling perturbation Ai a (l) scales as the bandwidth for a < whereas for a > 0, 
Fig. [2] the perturbation scales as A a /N a ^ 2 . 

triangles-curve in Figs. [2jand[3j where the PT-symmetry is exact. 
4. Chirality across the PT-symmetric Phase Boundary 

The time-evolution of a wave packet that is initially localized to a single site has 
been traditionally used to probe the degrees and signatures of PT-symmetry breaking 
in coupled optical waveguides [TBI HH]- 111 these systems, the wave function \ip(t)) 
denotes the single-transverse-mode electric field in each waveguide [27], the PT- 
symmetric impurities which absorb or emit the corresponding electromagnetic radiation 
are engineered, and the time evolution of an initially normalized wave packet is given by 
\ip{t)) = exp(—iH(X)t/h)\i/j(0)) where h is the scaled Planck's constant. The time- 
evolution operator exp(—iH(X)t/h) is not unitary irrespective of whether the VT- 
symmetry is exact or broken; therefore, in general, the net intensity I(t) = ^jl{j,t) 
is not one, where I(j,t) = \ (j\ip(t)) \ 2 denotes the site- and time-dependent intensity. 
This violation of unitarity is determined by the parity of number of lattice sites iV [26J , 
and the localized or extended nature of eigenstates [30]. In all cases, however, the 
net intensity increases monotonically across and exponentially past the PT-symmetric 
threshold [S]. 

The typical time-evolution of an initially localized state shows that apart from 
spreading across different sites, the wave packet undergoes a preferential clockwise 
or ant i- clockwise motion around the ring. We quantify this tendency, chirality, by a 
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dimensionless, Hermitian, discrete-momentum operator on a ring [22] . 

n(t) = m)\mt)) = -\ p^§^^ (5) 

where \ip{t)) = Y2j fj{t)\j) is a time-evolved wave function, and the normalization factor 
in the denominator is required due to its non-unitary time evolution. Note that due to 
the Cauchy-Schwartz inequality, the dimensionless momentum satisfies —1 < p^(t) < 1. 

When the Hamiltonian is Hermitian and the initial state is localized to a single site, 
the momentum p^(t) symmetrically oscillates about zero. As the impurity strength 
approaches the threshold value, 7 — > 7pt(^, «), Vipif) remains constant over long 
time intervals T ~ lOOiV ^> N where the time interval T is measured in units of 
2nh/A a [26]. We use the time-averaged momentum to denote this steady state value, 
Pipil) — Jo Pip(t')dt' /T, and choose the time interval T such that the steady-state value 
Pipi'y) is independent of it. 
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impurity strength Y/Y PT (m,a) impurity strength yly (m,a) 



Figure 4. Evolution of the dimensionless, average momentum (7) for N = 20 lattice 
with gain impurity +17 at position m — 2 (/i = 0.1), and initial wave packet locations 
too = 1 (blue circles dotted line) and too = 11 (red squares solid line). The left-hand 
panel corresponds to tunneling profile with a = 1, whereas the right-hand panel has 
results for a — 2. When 7 = the average momentum of the wave packet is zero. For 
small 7, hrst-order perturbation theory implies that ^(7) oc 7 with a slope that is 
dependent upon the initial state. In each case, the momentum shows a maximum at 
the threshold 7 = jpt and decreases monotonically on both sides of it [22], although 
the net intensity increases exponentially beyond the threshold [31 . 

Figure [4] shows the evolution of the dimensionless momentum ^(7) across the VT- 
symmetric phase boundary for an iV = 20 ring with tunneling functions a = 1 (left-hand 
panel) and a = 2 (right-hand panel). The initial positions of the wave packet are at 
m = 1 (blue circles dotted line) and m = 11 (red squares solid line). Note that when 
[i = m/N = 0.1, the threshold impurity strengths 7pt(/x)/A q , for a = 1 and a = 2 
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differ by an order of magnitude (Fig [2]). In all cases, the average momentum is zero 
when 7 = 0. Its linear increase with 7 at small j/jprpim, a) and the different slopes for 
different initial wave packet locations are both expected from a first-order perturbation 
theory. Remarkably, in all cases, p 1 p{'~f) reaches the same maximum possible value (of 
unity) at the PT-symmetric threshold and decreases monotonically beyond it even 
though the net intensity I(t) increases exponentially past the threshold. 

At this point, it is worthwhile to recall the corresponding results for a VT- 
symmetric ring with two (constant) tunnelings between the loss and gain impurities. 
In that case, the threshold strength 7pt(^) is independent of the impurity location m, 
the sign of the momentum ^(7) is determined by the path with the higher tunneling 
amplitude, and it reaches a universal maximum value of one at the threshold [26] . 




Figure 5. Dimcnsionlcss momentum results for an N = 22, a = 1 ring with initial 
wave packet at locations m = 1 (blue circles dotted line) and m,Q = 11 (red squares 
solid line). The left-hand panel has gain impurity at m = 8 (/z = 0.364) whereas the 
right-panel has gain impurity at location m = 7 (fi = 0.318). Due to the position- 
dependent tunneling Hamiltonian, analytical investigation of momentum dependence 
on gain impurity location m is not possible. These numerical results show that the 
sign and the maximum value of the momentum at the threshold are both dependent 
upon m, but not on the initial wave packet location mo. 

Figure [5] shows that for a ring with non-uniform, position-dependent tunneling 
profile t a , the behavior of the steady-state momentum is not as straightforward. Both 
panels present results for an iV = 22 ring with a = 1 and two initial wave packet 
locations, m = 1 (blue circles dotted line) and m = 11 (red squares solid line). The 
left-hand panel has the gain impurity at location m = 8 (fi = 0.364) and it shows that 
the sign of the steady-state momentum is now negative. The right-hand panel has gain 
impurity at location m = 7 (fi = 0.318) and it shows that the maximum value attained 
by the steady-state momentum is not unity. These results imply that, contrary to the 
two-tunneling ring, sign of the momentum and the maximum value it attains at the 
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threshold are both dependent upon the impurity location. We emphasize, however, that 
these results are independent of the initial wave packet location, and the qualitative 
behavior ofp^j) across the VT -symmetric threshold is identical in all cases. 

5. Discussion 

In this paper, we have numerically explored the "PT-symmetric phase diagram for 
a finite, nonuniform PT-symmetric lattice as a function of its boundary conditions. 
We have shown that for experimentally relevant lattice sizes, the differences between 
properties of an open chain and a ring are nontrivial, particularly for tunneling profiles 
t a (k) with a < 0. Generically, we found that for a > 0, the PT-symmetric phase 
in a ring is weaker than its counterpart in an open chain. In contrast, when a < 
the PT-symmetric phase in a ring is substantially strengthened. Since for a < the 
tunneling perturbation that is required to change an open chain into a ring is comparable 
to the bandwidth of the open chain, the strengthening of the PT-symmetric phase is 
reasonable. 

We have shown that the VT symmetry breaking is accompanied by a qualitatively 
universal behavior of dimensionless, average momentum ^(7): the momentum is zero 
when 7 = 0, increases linearly with 7, and its magnitude reaches maximum at the PT- 
breaking threshold, accompanied by monotonic decay on both sides of the threshold. 
We have also found that, in contrast with the two-tunneling model [22], here, the sign 
of the momentum and its maximum value at the threshold are dependent upon the 
impurity locations, but not on the initial wave packet location. 

These numerical results raise several questions. For an open chain, the location of 
the pair of eigenvalues that become degenerate and then complex is uniquely determined 
by the impurity location fi irrespective of the value of a [26J; no such claim seems 
possible for the VT symmetry breaking in a ring. The universal presence of the peak in 
|p^(7)| at the threshold 7 = jpri^, ot) suggests that it may be driven by the exceptional 
point at the threshold where two eigenvalues become degenerate and the corresponding 
eigenvectors become parallel [32J. However, lacking analytical methods for the non- 
uniform tunneling ring, a systematic numerical investigation of the dependence of the 
sign and the maximum value of the chirality remains an open problem. 
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